Conserving and gapless approximations for the composite bosons in terms of the 

constituent fermions 
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A long-standing problem with the many-body approximations for interacting condensed bosons 
has been the dichotomy between the "conserving" and "gapless" approximations, which either obey 
the conservations laws or satisfy the Hugenholtz-Pines condition for a gapless excitation spectrum, 
in the order. It is here shown that such a dichotomy does not exist for a system of composite bosons, 
which form as bound-fermion pairs in the strong-coupling limit of the fermionic attraction. By 
starting from the constituent fermions, for which conserving approximations can be constructed for 
any value of the mutual attraction according to the Baym-Kadanoff prescriptions, it is shown that 
these approximations also result in a gapless excitation spectrum for the boson-like propagators 
in the broken-symmetry phase. This holds provided the corresponding equations for the fermionic 
single- and two-particle Green's functions are solved self-consistently. 

PACS numbers: PACS numbers: 03.75.Ss, 03.75.Hh, 05.30.Jp 



Many-body decriptions of a system of interacting con- 
sensed bosons have long been known [ij to fall into either 
one of two classes of approximation schemes, which are 
alternatively consistent with the conservation laws (con- 
serving approximations) or with the absence of a gap 
in the elementary excitations spectrum (gapless approx- 
imations). Having to choose between these two types of 
approximations constitutes a shortcoming of the many- 
body theory for condensed bosons, as one would rather 
like to deal with approximations which are conserving 
and gapless at the same time. 

These approximation schemes have been conceived for 
interacting-boson systems like helium, for which the in- 
ternal fermionic structure is immaterial, due to the large 
amount of energy required to excite the internal fermionic 
degrees of freedom compared with the energy scales of the 
experiments. 

Recent experimental advances with ultracold trapped 
Fermi atoms have made it possible to produce systems 
of composite bosons (dimers), whose binding energy is 
comparable with the energy and temperature involved in 
the experiments |2|- The Bose- Einstein condensation of 
the dimers has also been detected |3|- In these systems, 
the internal fermionic structure is definitely relevant, as 
the binding energy of the dimers can be tuned across 
threshold via a Fano-Feshbach resonance |j|, sweeping 
from bound to unbound fermions and viceversa. 

For these systems, it appears appropriate to construct 
the dynamical propagators of the condensed composite 
bosons in terms of the constituent fermions, by follow- 
ing the progressive quenching of the fermionic degrees 
of freedom as the fermionic attraction is increased. For 
the constituent fermions, it is known that conserving ap- 
proximations can be constructed for any value of the mu- 
tual attraction even in the broken-symmetry (superfluid) 
phase, via the Baym-Kadanoff prescriptions ^5;, ;6J which 
require the self-consistent solution of the equations for 



the single- and two-particle Green's functions. In this 
way, conservation laws can be regarded to be fulfilled 
not only in terms of the constituent fermions but also 
in ternrs of the composite bosons, when the fermionic 
attraction gets sufficiently strong that the fermionic de- 
grees of freedom (internal to the composite bosons) are 
progressively quenched. 

The question is whether such conserving approxinra- 
tions (for the constituent fermions and therefore for the 
composite bosons) could also result into a gapless ex- 
citation spectrum for the propagators of the composite 
bosons in the broken-symmetry phase, the gapless condi- 
tion being required ongeneral grounds by the occurrence 
of a Goldstone mode [3| . [The case of a homogeneous sys- 
tem will be specifically considered when discussing the 
absence of a gap in the bosonic excitation spectrum [g .] 

Purpose of this paper is to show that a given fermionic 
conserving approximation also results in a gapless ex- 
citation spectrum for the boson- like propagators. The 
long-standing dichotomy between conserving and gapless 
approxinrations can thus be apparently overcome, when 
the bosons themselves are treated at a more fundamental 
level in terms of the constituent fermions. 

We begin by generalizing to the composite bosons the 
theorem of Hugenholtz-Pines for ordinary bosons, by fol- 
lowing the treatment of Ref . llj in terms of formally exact 
propagators (see also Ref. 0). Fermionic conserving ap- 
proximations to select approximate propagators for the 
composite bosons will be considered later. To this end, 
we define a bosonic-like field operator 

^B{r) - ldpcb{p)ij^{r - p/2)V'T(r + p/2) (1) 

for any value of the fermionic coupling, where -tpa- (r) is a 
fermionic field operator with spin a. When the fermionic 
attraction is sufficiently strong, on physical grounds the 
(real and normalized) function (l){p) can be taken as the 



bound solution of the associated two-body problem. At 
weaker coupling, a precise choice of (j){p) is not required. 
For instance, it could be taken as the solution of a gen- 
eralized Cooper problem, whereby the Fermi energy is 
replaced by the coupling- (and temperature-) dependent 
fermionic chemical potential ^. To break the gauge sym- 
metry, the bosonic order parameter Q;(r) = (\I/B(r)),, 
is defined as the thermal average of the operator (Q) 
within the restricted (77) ensemble of Ref. |lj. With the 
Nambu representation for the fermionic field operator 
(^i(r) = ■(/'t('')' ^2(1") = "^Ki")); this thermal average can 
then be expressed in terms of the anomalous fermionic 
single-particle Green's function Q12 

(*^(r-f )*i(r+f )), = ei2(r+f ,r-|;r = 0") (2) 

with imaginary time r. 

In what follows, it is convenient to consider a general- 
ized fermionic single-particle Green's function: 



5(1,1') 



{Tr[S]) 



(3) 



with the notation 1 = (ri, ri, £1) in terms of the Nambu 
spinor component £. Here, Tr is the imaginary-time- 
ordering operator, (■ • ■) is a thermal average taken with 
the system grand-canonical Hamiltonian K = H — fiN, 
\I'(1) = exp{7^Ti}*£j(ri) exp{— Xti}, and the operator 
S = exp{- / dll'^^l)U{l, l')^(l')} contains the source 
term 



L/(l,l') 



C/,(ri,ri-;Ti) -C/„(ri, n-; n) J ^^^^ ^^' > 

(4) 
with a normal (Un) and a superfluid (Us) component. 
[The normal component C/„ will not be needed in the 
following, while the superfluid component Us will be al- 
lowed to vanish at the end of the calculation.] In the 
static case (when U does not depend on the imaginary 
time), the generalized definition (|2Jl concides with the 
ordinary definition [like Eq.Q] within the 77— ensemble. 
With the definition Q, the bosonic order parameter is 
generalized as follows: 



a(r) 



dp(t>{p)gi2(r- 



P 

2' 



(5) 



Suppose now that [/s(r,r';T) is varied by a small 
uniform change of phase (5$, such that (5C/s(r, r'; r) ^ 
«(5<I>[/s(r, r'; r). This change can be reabsorbed by a 
canonical tranformation of the fermionic field operators, 
so that the corresponding change of the order parameter 
© is given by ^a(r) = —iS^air) to the leading order in 
(5$. The change (5a(r) can be calculated alternatively via 
the definitions (jSJ and Q, by performing the functional 
derivative of (j2Jl with respect to a variation of Us- One 
obtains: 



5a{v) = dp(j>{p) / dr2 dr'2 / dT2 



L{v + |, r, 1; r2, T2, 2; r - f , t+, 2; r'^,T+ , 1) 

y.5Us{V2,Y'^]T2)* 



-,r, I;r2,r2,l;r 
x,5f/,(r'2,r2;T2)] 



-L(r + |, r, 1; ra, T2, 1; r - |, t+, 2; r^, t+ ,2) 



(6) 



Here, i(l,2,l',2') = ^2(1, 2, 1', 2') - ^(1, l')a(2, 2') is 
the two-particle correlation function expressed in terms 
of the generalized fermionic two-particle Green's function 



a2(l,2,l',2') 



(r,[5^(l)^(2)^t(2')^t(i/)]) 
{Tr[S]) 



, (7) 



and is obtained via the functional derivative 
L(l,2,l',2') = -<5g(l,l')/^C/(2',2). 

By a similar token, for the adjoint a(r)* of a(y) one 
obtains (5a (r)* = i(5<I'a(r)*, as well as 

(5a(r)* = dp(f>{p) / dr2 dr2 / dT2 

L(r-|,r,2;r2,T2,2;r+|,r+,l;r^,T2+,l) 
x6Us{r2,r2;T2y 



P P 

|,r,2;r2,r2,l;r ' ^ 

xSUs{r'2,r2;T2)] 



+ L{v - |, r, 2; ra, T2, 1; r + |, t+, 1; r^, T2+, 2) 



(8) 



in the place of jnj. [The quantity a(r)* is defined like 
in Eq.(|SJl with ^21 replacing Q12, and coincides with the 
complex conjugate of a(r) in the static case.] 

The static and uniform limit of the above results can 
be considered at this point. Accordingly, without loss 
of generality we let C/s(r,r';T) — > C/s0(|r — r'|) in both 
Eqs. jnj and ((HJ, where Us is a (complex) constant and 
(/)(r) the same function of Eq. (QJ . We also introduce the 
Fourier representation: 



(9) 



dq 1 



gi(p+q)-ri gjp'-r2 g-ip-r'i g-i(p'+q)-r2 



(27r)3/3 2-' 

Here, p, p', and q are wave vectors, a;„ = (2n + l)7r//3 
(n integer) is a fermionic Matsubara frequency (/3 be- 
ing the inverse temperature), Qi, = i'Kv/P [v inte- 
ger) a bosonic Matsubara frequency, and p — (p,Li;„), 
p' — (p'^ujni), and q = (q, ill/) is a four-vector notation. 
By straightforward manipulations of Eqs. © and JHJ, 
and by recalling the identities 5a{Y) = — i(5(f>a(r) and 
6Us{v,v'\t) = ii5$[/s(r,r'; r) (plus their adjoints), one 
ends up with the matrix equation 



G\\{q^Q) G\l{q 
Glliq^O) Glliq 



0) 
0) 



Us 
-Us 



(10) 



for U and a real, with the notation 



GZliq) 






^i4 



(p + q/2)(/.(p' + q/2)4^;^(p,p'; g) . (11) 



It can be readily verified that, in the limit Us — > 0, the 
definition 1)11(1 corresponds to the four possible (normal 
and anomalous) bosonic-like propagators which can be 
constructed with the operator ^ and its adjoint. 

Before letting C/j — *■ in Ea. (|10l) . it is convenient to 
introduce the inverse of the matrix on its right-hand side 
and write 

GM(g^O) GlUq^O)\ _ 1 f D -B 

Glliq ^ 0) Glliq ^ 0) J AD-BC\-C A 

(12) 
Matrix inversion of Eq. (|10|l then yields the conditions 

A-B = , C-D = , (13) 

in order to have a finite value for the order parameter 
a in the limit Us — > 0. These conditions are not inde- 
pendent from each other, since one can prove on general 
ground from time-reversal invariance that A — D and 
B ~ C. The denominator in Ea. H12() then reduces to 
AD - BC = {A - B){A + B), and vanishes owing to 
(|13|l . This implies that the bosonic-like propagators Hll|) 
are singular when g — > 0, irrespective of the value of the 
fermionic coupling. In the present context, the condition 
A — B — corresponds to the Hugenholtz-Pines theorem 
for ordinary bosons 9]. It implies, in particular, that the 
composite bosons, which form when the fermionic attrac- 
tion is strong enough, have a gapless spectrum. 

All considerations made so far hold for the ex- 
act fermionic single- [Eq.||2J)] and two-particle [Eq.(I7|)] 
Green's functions. The crucial point to derive the re- 
sult Hl()|l was that the single- and two-particle fermionic 
Green's functions are related to each other via a func- 
tional differentiation in the presence of the external po- 
tential U. When dealing with fermions, it is standard 
practice to explore this relation further [a, |6j , by exploit- 
ing the Dyson's equation with the self-energy S: 

- g-\l, 2) = ^5(1, 2) + M(l, 2) + [/(I, 2) + E(l, 2) . 

(14) 
In this expression, M(l,2) — rf j, (/i(ri) — ^)6{rx — 
r2) '^(■''i — '''2) where r^ is a Pauli matrix and /i(r) is the 
single-particle Hamiltonian (which, in general, includes 
an external static potential). The two-particle correla- 
tion function L is correspondingly obtained as: 



-L(l,2,l',2' 
--- fd34g{l,3) 



■/ 



d3456g(l, 3)^(6,1') 



SU{2',2) 

5g-H3A) 

dU{2',2) 

.,5S(3,6) 



a(4,l')=e(l,2')e(2,l') 
(-)L(4,2,5,2').(15) 



It thus satisfies the Bethe-Salpeter equation, with kernel 
6I]/5Q related to the kernel E of the Dyson's equation 
(tn)l . The limit U -^ can be taken in Eqs. ifT^ and 
(|15|l whenever appropriate. 

Selection of an approximate fermionic many-body the- 
ory starts with an approximate choice of the func- 
tional form of the self-energy E in terms of Q (and 
of the two-body interaction). The equations H14|l and 
()15|l are then solved self-consistently, with the respec- 
tive approximate kernels E and dT,/dQ. In addition, 
Eq.l^ni) implies that Eq.® (and its adjoint ((SJ)) holds 
even for the approximate theory, since L (with the ap- 
proximate kernel 5Y,/dQ) still represents the functional 
derivative of Q with respect to U. The alternative re- 
sult 5a{r) — — i(5$a(r) is instead obtained in the ap- 
proximate theory by noting that, under the transforma- 
tion t/s(r, r';T) -^ e*''*[/s(r,r'; r), the approximate off- 
diagonal single-particle Green's function Q12 of Eq.jSJ 



trasforms as Q12 



i5<S> 



Q12. As a consequence, the result 



<5e(4,5) 



()13|l follows even within the approximate theory, imply- 
ing a gapless spectrum. 

Conserving approximations for fermions are similarly 
based on Eqs. lfn|l and ((T3|l . which are solved self- 
consistently for a given approximate choice of E. In this 
case, the self-energy has to be chosen appropriately, to 
comply with the requirements of local number conserva- 
tion and gauge invariance jlflj. It is then required that 
the symmetry property L(l,2,l',2') = L(2,l,2',l') is 
satisfied by the approximate L. To this end, it is suf- 
ficient that the approximate kernel 5T,/5Q of Eq. 1(1 5|l 
satisfies the same symmetry property. This property is, 
in turn, met by any $— derivable approximation for the 
self-energy E of Eq.(ini), whereby E(l,2) = 5<^/5g{2, 1) 
is obtained from an approximate functional <&. j^ Atten- 
tion must be paid to the fact that a given choice for E 
may not meet this requirement, unless certain diagrams 
for E are taken together j^ |ll| . 

It is now evident that a conserving approximation for 
the constituent fermions, which holds for any value of 
their mutual attraction, will also result in a conserv- 
ing approximation for the composite bosons that form 
in strong coupling. The same fermionic conserving ap- 
proximation will further result in a gapless spectrum 
for the composite bosons, as the same requirement for 
Eqs. 1(14(1 and 1(15(1 to be simultaneously self-consistently 
satisfied applies to both (conserving and gapless) proce- 
dures. This proves our claim. From the above considera- 
tions, it is also clear that the requirements for a fermionic 
approximation to be conserving at any given coupling are 
more stringent than the absence of a gap in the excitation 
spectrum of the composite bosons in strong coupling. 

A well-known example of a fermionic conserving ap- 
proximation, which results for any coupling in a gapless 
spectrum for the collective mode associated with the bro- 
ken symmetry [10|, is the ECS (off-diagonal) approxima- 
tion for E, shown in Fig. 1(a) with the associated po- 
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FIG. 1: Self-energy E derived from the potential $ with 
the associated kernel ST,/SQ, for (a) the BCS approximation 
and (b) the t-matrix approximation in the broken-symmetry 
phase. Full lines represent fermionic (self-consistent) single- 
particle Green's functions, with the arrows pointing from 
the second to the first argument; broken lines represent the 
fermionic interaction potential. 



similarly constructed by suitable choices of the fermionic 
self-energy. 

In conclusion, we have shown that a given conserv- 
ing approximation for the constituent fermions also re- 
sults into a gapless spectrum for the composite bosons. 
By following the formation of the bosons from the con- 
stituent fermions as the fermionic attraction is progres- 
sively increased, a long-standing (conserving vs gapless) 
dichotomy can thus be resolved, at least at a formal 
level. Although the self-consistent solution of the equa- 
tions determining the fermionic single- and two-particle 
Green's functions might, in general, involve considerable 
numerical labor, enforcing the fermionic conserving crite- 
ria proves per se sufficient to get a gapless bosonic spec- 
trum. 

We are indebted to A. Fetter, F. lachello, and G. 
Morandi for a critical reading of the manuscript. This 
work was supported by the Italian MIUR (contract Cofin- 
2003 "Complex Systems and Many-Body Problems"). 



tential <& and kernel 5T,/5Q of the Bethe-Salpeter equa- 
tion H15|) . In this case, the self-consistent solution of the 
Dyson's equation ifT^ reduces to the solution of the BGS 
gap equation; in addition, this equation coincides with 
the condition A — B ~Q oi Eq. H13|) which guarantees the 
absence of a gap in the bosonic excitation spectrum. In 
this simple case, therefore, the self-consistent solution of 
the BCS gap equation is sufficient for the approximation 
to be conserving and gapless. More generally, separate 
solutions of the equations H14|) and (|15|l are required for 
the approximation to be conserving and gapless. For in- 
stance, to the self-energy E within the fermionic t-matrix 
approximation in the broken-symmetry phase |l2l | , shown 
in Fig. 1(b) with the associated potential $, there corre- 
spond three distinct contributions to the kernel 5T,/5Q, 
also shown in Fig. 1(b). When considering the BCS and 
t-matrix approximations for S together, to get a gapless 
spectrum it is thus not enough to solve self-consistently 
the Dyson's equation for Q with both self-energy con- 
tributions, if one solves at the same time the Bethe- 
Salpeter equation with only the BCS contribution to 
kernel 5Yi/5Q. By doing so, one would, in fact, omit 
the three contributions to 5Y,/5Q depicted in Fig. 1(b), 
whose presence is required by conserving criteria. Ad- 
ditional conserving and gapless approximations can be 
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